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QCD can be formulated using any gauge group. One particular interesting choice is to replace 
SU(3) by the exceptional group G2. Conceptually, this group is the simplest group with a trivial 
center. It thus permits to study the conjectured relevance of center degrees of freedom for QCD. 
Practically, since all its representation are real, it is possible to perform lattice simulations for 
this theory also at finite baryon densities. It is thus an excellent environment to test methods 
and to investigate general properties of gauge theories at finite densities. We review the status 
of our understanding of gauge theories with the gauge group G2, including Yang-Mills theory, 
Yang-Mills-Higgs theory, and QCD both in the vacuum and in the phase diagram. 
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1. Introduction 



Gauge theories can in principle, up to anomalies, be formulated for all simple Lie groups. This 
property has been used often to gain insight into structures or simplify calculations. One salient 
example is the large-A^ limit in QCD. Another option is to use the exceptional group G2, leading to 
G 2 QCD. 

The proposal to make this replacement was made [1] to understand the role of the center of 
the gauge group, which was long assumed to play a central role for many of the salient features of 
QCD, especially confinement. However, the detailed investigations, to be presented in section ||, 
showed that most of these features are also present in the G2 case. 

Besides these conceptual questions concerning the center, another property of G2 QCD is 
interesting from a practical point of view. Since all its representations are real, no sign problem 
arises when simulating G2 QCD with dynamical fermions. It is thus possible to investigate the 
whole phase diagram of the theory using lattice calculations [2]. G2 QCD is so far the theory most 
similar to QCD where this is possible in the continuum limit. The resulting phase diagram [2] is 
rather similar to the one obtained in other such theories, like QCD with gauge group SU(2) (QC2D) 
[3-6] or QCD in the strong coupling limit [7, 8]. Thus, G2 QCD offers another perspective on the 
QCD phase diagram. This will be detailed in sections |] and |5[ 

It is, of course, an interesting question whether there can be established any direct connection 
between the G2 case and the SU(3) world. Breaking the G2 gauge group using a Higgs field works 
for the Yang-Mills case [1], as briefly outlined in section |3|, but it is yet not clear whether this is 
also possible in the QCD case. 

Thus, gauge theories with gauge group G2 are very interesting from many perspectives, as will 
be summarized in section ||. However, most investigations are yet on a qualitative and exploratory 
level, and many interesting questions have not even been addressed yet. 



2. Yang-Mills theory 
2.1 Zero temperature 

The simplest realization of a gauge theory with the gauge group G2 is Yang-Mills theory. 
Since the adjoint representation of G2 is 14-dimensional, there are 14 gluons. Using the Macfarlane 
representation [9] a G2 link (or group element) U in the 7-dimensional fundamental representation 
can be written as 

(u 0' 
U=Z 10 
^0 u* 

where Z is a 7-dimensional representation of S 6 and u is an element of SU(3). Thus, 8 of the gluons 
can be considered loosely as 'SU(3)'-like. This will become important in section ||. Due to this 
explicit SU(3) subgroup, lattice simulations of a G2 theory are straightforward but expensive, see 
[2, 10-12] for the algorithms employed here. 

G2 is the smallest rank 2 gauge group with a trivial center. As a consequence, every funda- 
mental charge can be screened by three adjoint charges, and thus there is no infinitely rising Wilson 
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Figure 1: The Wilson potential V divide by the scale fi for different representations ffl (left) and its string 
breaking, compared to hybrid masses for two representations (right), both in three dimensions, from [12]. 
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Figure 2: The minimal Landau-gauge gluon propagator D (left panel) and running coupling a (right panel) 
of G2 Yang-Mills theory compared to SU(3) Yang-Mills theory in three dimensions as a function of momen- 
tum p, from [2]. 



potential, and thus no confinement in the sense of a Wilson area law [1]. However, in practice the 
corresponding Polyakov loop is found to be very small at zero temperature, and in fact only upper 
bounds are known, though it follows that it must be non-zero. In fact, at intermediate distances a 
linear rising Wilson potential [13, 14], including a characteristic Casimir scaling [12, 15], is found. 
Thus, a string appears in the same way as in QCD with dynamical quarks, up to a distance where 
string -breaking sets in [12]. Hence, G2 Yang -Mills theory is in the same sense (non-)confining as 
is QCD. These facts are illustrated in figure |]. Of course, since the theory has no anomaly, it is still 
a well-defined theory, with only colorless asymptotic states [1, 13], like glueballs [12, 16]. 

It is thus an interesting question what the effective degrees of freedom are. On the level of the 
elementary particles, the gluons, no qualitative, and little quantitative difference is found [10, 17]. 
This also manifests itself in a qualitative similar running coupling, even in the far infrared, see 
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figure ^ . Thus, at the level of gluons, there is no distinct difference. 

Another set of effective degrees of freedom often used in Yang-Mills theory are topological 
ones. Similarly, for G2 Yang-Mills theory vortices [14], monopoles [19], dyons [20], and instantons 
[21] have been constructed. Using lattice simulations and cooling, it is indeed possible to verify 
the existence of topological lumps, which are associated with action lumps and a non-vanishing 
topological susceptibility of roughly (150 MeV) [21], though yet with large systematic errors. 
Though there exist differences in details, e. g. vortices are not associated with a center [14], the 
salient features of these topological excitations are close to the ones in ordinary SU(A0 Yang-Mills 
theory. As one can expect from these observations, chiral symmetry is broken in the vacuum in the 
same way as in ordinary Yang-Mills theory [18]. 

Thus in total, G2 Yang-Mills theory in the vacuum is very similar to S\J(N) Yang-Mills theo- 
ries. 



2.2 Finite temperature 

Since the finite-temperature phase uansition in SU(A0 Yang-Mills theories is associated with a 
center-symmetry breaking/restoring phase transition, it was originally anticipated [1] that there will 
not be a phase transition in G2 Yang-Mills theory, though the gluonic sector suggested otherwise 
[22]. Lattice simulations then indeed found a strong first-order phase transition in G2 Yang-Mills 
theory [13, 14, 23] using the free energy. However, in practice this is non-trivial due to a bulk 
transition requiring rather fine lattices [13, 23]. This phase transition is also reflected in the behavior 
of glueballs [16]. 

Amazingly, though not being an order parameter, the Polyakov loop also reflected this phase 
transition. In fact, it is possible to use the Polyakov loops in various representations to describe the 
phase structure of G2 Yang-Mills theory rather accurately [24]. One of the main reasons seems to 
be that though there is no genuine center symmetry, a distorted three-fold structure is still preserved 
by G2, which, when breaking the theory down to SU(3), yields the center symmetry, see section [| 
below. 

This alone is already in remarkable agreement to ordinary Yang-Mills theory. But the similar- 
ities are even more pronounced. Since all representations are real, it would have been possible that 
the chiral transition, as is the case for the adjoint chiral condensate in S\J(N) [25, 26], would not 
show a phase transition or only at a much higher transition temperature. This is not the case, and, 
within lattice resolution, the chiral condensate shows a response precisely at the same temperature 
as the Polyakov loop and the free energy [18]. As would be naively expected from the comparison 
to SU(A0 Yang-Mills theory, it is then also found that the topological properties change at the phase 
transition [21], especially the topological susceptibility drops. 

The resulting phase diagram is shown in figure |3[ The first order nature is visible, though it 
requires a detailed study of scaling properties to ascertain it [23]. Thus, from the point of view of 
the phase diagram G2 Yang-Mills theory behaves very similar to the SU(3) case, even though the 
phase transition is not related to a symmetry. This is one of the reasons why Yang-Mills theory is 



well suited as a stand-in for QCD thermodynamics, as discussed in section 4.2. The reason for the 



'For all results for Yang-Mills theory, the scale has been set by giving the intermediate distance fundamental string- 
tension a value of (440 MeV) 2 [10, 18]. 
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Figure 3: The phase diagram of G2 Yang-Mills theory. The critical temperature is taken from [23], the 
Polyakov loop and chiral condensate from [18] and the topological susceptibility from [21]. 

existence of this similarity is besides the approximate three-fold structure [24] the fact that the size 
of the gauge group appears to be more relevant for the phase structure than the center of the group 
[13,27]. 



3. Yang-Mills-Higgs theory 

One of the interesting features of G2 is that it has SU(3) as a sub-group. Thus, it appears 
possible to somehow hide the S 6 part of the gauge group using the Higgs mechanism such that 
just SU(3) remains. In fact, it turns out that a single fundamental Higgs field is sufficient for this 
purpose [1, 11]. In such a more complicated theory it is possible to follow the phase structure 
at finite temperature, and map a phase diagram in the temperature-Higgs mass plane at infinite 
four-Higgs coupling [11], as shown in figure^. 

The phase structure is rather intricate at intermediate values of the couplings. Given the large 
systematic uncertainties encountered in such theories [28] a definite answer will remain hard to 
find. However, this question is highly relevant: If a continuous connection between the SU(3)-like 
domain and the G2-like domain exists this would have significant implications for the physics of 
both theories. 

The situation becomes much more complicated when introducing (fermionic) matter fields 
into the theory [1]. In this case, a hiding with just one Higgs field will inevitably lead to an SU(3) 
theory with the matter fields in the wrong representation, in particular to real matter fields. Since 
the natural question is, whether a connection to ordinary QCD is possible, the hiding or breaking 
mechanism must complexify the matter fields to lead to the inequivalent fundamental and anti- 
fundamental representations of QCD. This will likely only be possible, if at all, by manipulating 
the theory on the level of Weyl fermions, a topic under current investigation [29]. 
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Figure 4: The phase diagram of G2 Yang-Mills-Higgs theory, from [11], as a function of gauge coupling 
and Higgs hopping parameter at finite temperature. 



4. G 2 QCD 

4.1 Vacuum structure 

When adding Nf fundamental fermions to G2 Yang-Mills theory one arrives at G2 QCD. The 
vacuum structure of this theory is yet little explored [1, 2], but has a number of highly interesting 
features. The first concerns the spectrum. Due to the group structure, there exists a richer set of 
color-neutral bound states than in QCD [1], both of fermionic and bosonic type. In the boson sector 
there are as in QCD the glueballs and mesons. In addition, there are also diquarks, since due to the 
reality of the G2 representations such states are color-neutral, different from QCD, but similar to 
QC2D. In addition, there are also tetraquarks and heptaquarks consisting out of four and six quarks. 
Besides these bosonic hadrons there are also fermionic ones. Most notably the hybrid, consisting 
out of one quark and three gluons, but also a nucleon from three quarks, as well as pentaquarks and 
heptaquarks from five and seven quarks. 

The mass hierarchy of these states will depend strongly on the masses of the quarks, even 
for degenerate flavors. E. g., at heavy quark mass the hybrid will be the lightest particle in the 
fermionic sector, while the nucleon is expected to take over this role at low quark masses, but will 
still be heavier than the diquark or mesons. The details of these hierarchy are a dynamical problem. 
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Figure 5: The connected part of the diquark/scalar meson correlator with a mass fit (left panel) and the 
masses for the diquark and the pion as a function of the gauge coupling (right panel), from [2]. The lattice 
spacing is strongly-dependent on the lattice parameters. 



These bound states are also influenced by the pattern of chiral symmetry breaking. Due to the 
reality of quarks, G2 has, similarly to QC2D, an enlarged chiral symmetry of U(2Nf) [1, 2, 30, 31]. 
This symmetry can be viewed as a flavor symmetry on the level of the Weyl fermions. Of this 
symmetry an axial U(l) is expected to be broken in the same way as in ordinary QCD by the axial 
anomaly. Taking for the following a single flavor leaves, in contrast to QCD, still an SU(2) chiral 
symmetry. This symmetry is spontaneously broken [2], like in the quenched case [18], leaving only 
an U(l) intact. This conserved U(l) can then be associated with a baryon number. The Goldstone 
bosons of this breaking are then expected to be two diquarks [2], just like in QC2D [5]. These two 
diquarks represent a flavor-doublet on the level of Weyl fermions. 

In numerical simulations this is rater hard to identify, as for one Dirac flavor the scalar and 
the diquarks only differ by disconnected contributions. Furthermore, it turns out that G2 QCD in 
the range of accessible parameters is very sensitive to both the gauge coupling and the hopping 
parameter, and has at rather low gauge coupling at fixed hopping parameters already a transition 
into an unphysical phase [29], possibly an Aoki-like phase. Nonetheless, mass determinations are 
possible, as is demonstrated for Nf = 1 in figure |5[ The determination of the vacuum spectrum is 
thus a challenging task, even at a qualitative level, and an ongoing project [29]. Especially the mass 
of the nucleon is relevant, when one turns to the phase diagram. 

4.2 Phase diagram 

Due to reality of the representations and the enlarged chiral symmetry, the whole phase dia- 
gram for the Nf = 1 case is both accessible in lattice simulations and relevant. Even besides the fact 
that G2 QCD is an interesting theory on an intellectual level, there is a number of features which 
makes it also highly relevant on the level of applications in the continuum limit. First of all, as de- 
scribed in section |2[ the theory is in the quenched limit very similar to SU(A0 Yang-Mills theories, 
in contrast to theories with adjoint matter [25, 26, 31]. Furthermore, the theory has nucleons, and 
in general fermionic baryons, and thus also nuclei. Hadronic Pauli effects at intermediate densities 
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Figure 6: The G2 QCD phase diagram for one flavor of quarks. The left panel shows the (unrenormalized) 
Polyakov loop, the middle panel the normalized chiral condensate, and the right panel the Baryon density, 
normalized to the saturation density of 14 quarks/lattice site. For details and simulation parameters, see [2]. 
Note that at /XQ uar k ~ 1 GeV the system starts to become dominated by systematic effects [2]. 

will thus play a role, in contrast to QC2D [3-6]. No other gauge theory with this combination of 
features has yet been simulated on a lattice, except without continuum limit [7, 8]. 

This provides the possibility of a number of unprecedented tests of lattice approaches to finite 
density QCD. It is possible to test explicitly to which extent investigations using analytical con- 
tinuation in imaginary or isospin chemical potential work (see e. g. [32, 33]), and whether Taylor 
expansions (see e. g. [34, 35]), Lee- Yang zeros (see e. g. [36]), or other methods (see e. g. [37]) are 
reliable tools. 

Furthermore, and possibly even more important, the G2 QCD lattice phase diagram provides 
new benchmarks for both models [38-40] and continuum methods, in the latter case especially 
functional methods [40-42]. Furthermore, if breaking G2 QCD to ordinary QCD should be possi- 
ble, this would be even more helpful, though, of course, at some point the sign problem will prevent 
a simulation of QCD. 

The first step in this program is provided by a proof-of -principle showing the accessibility of 
the phase diagram in lattice calculations, see figure |p [2]. Though so far at a qualitative level, 
it shows already a structure close to the expected one, including indications [29] of a silver-blaze 
point [43], see also section ||| A more quantitative description will require more detailed calcu- 
lations, in particular concerning systematic errors [29]. Nonetheless, the theory shows the low- 

2 The scale is here chosen to get a zero-density transition temperature of about 160 MeV, and the first excited meson 
state is used to set the scale scale. This procedure [2] is strongly affected by systematic errors, and will be improved 
[29]. 



8 



G2 gauge theories 



Bjorn H. Wellegehausen 




Figure 7: The baryon density (left panel) of G2 QCD, compared to that of continuum and lattice [3] Stefan- 
Boltzmann results with the same mass or for massless quarks, and to leading order chiral perturbation theory 
[3, 30] with coefficients fitted to the G2 case at intermediate densities. The middle panel shows the cor- 
responding ratios (note the logarithmic scale), and the right-hand panel the integrated equation of state, 
normalized to the continuum Stefan-Boltzmann case. All results unpublished from [29]. The value of the 
lattice constant a is approximately 0.2 fm. 



temperature, low-density ordinary phase, has a transition, likely a cross-over, to a high-temperature 
phase, and also a transition at finite density. Whether any of these are phase transitions remains 
to be seen, but so far the finite-density transitions are stronger. Also, first signals of additional 
structure at zero temperature have been observed [29], which may correspond to various phase 
also observed in QC2D [5, 6]. However, more details studies, especially of systematic effects are 
necessary before definite statements can be made. 

Finally, an interesting question is, to which extent the low-temperature case is simple, so that e. 
g., quasi-particle models would be a good description. For this purpose, a comparison to a system 
of free quarks and to chiral perturbation theory is shown in figure 0. While the high-density region, 
which is dominated by lattice artifacts [2], is rather well described by the corresponding free lattice 
system of quarks, this is not the case at low densities. Here, the equation of state is much more 
similar to lattice or continuum versions of a gas of free massless quarks, instead of massive ones, 
though the deviations are still very large at the smallest densities. At the same time, at least leading- 
order chiral perturbation theory is not able to reproduce even qualitatively the physics of G2 QCD. 
Thus, non-trivial effects play a dominant role at densities below a\i ss 0.5, which translates in this 
case to roughly 500 MeV of quark chemical potential. In this region, highly non-trivial effects have 
to be dealt with. 



5. Results on a smaller lattice 

Since many of the investigations above are limited by the number of different lattice settings 
which can be simulated, the use of smaller lattices may help in mapping the phase diagram on a finer 
grid. However, due to the unphysical bulk transition it is not possible to study the full phase diagram 
on smaller lattices, especially at finite temperature on lattices with N t < 5. Nevertheless, at zero 
temperature G2 QCD is investigated on a 8 X 16 lattice in the parameter region j8 = 0.90. . . 1.10 
and jU = ... 2. The monopole density is already sufficiently small, such that the system stays 
outside the bulk phase for all values of j8 and jj.. In Fig. |8] the diquark and nucleon (proton) 



9 



G2 gauge theories 



Bjorn H. Wellegehausen 



4.0 
3.5 
3.0 
2.5 
m 2.0 
1.5 
1.0 
0.5 
0. 



~i 1 1 r 

m diqii<irk 1 * I 



1 * X 



(.90 0.95 



_I L_ 



1.00 

p 



1.05 1.10 




0.15 



0.90 0.95 1.00 1.05 1.10 



Figure 8: Diquark and nucleon mass (left panel) and its ratio (right panel) on a 8 x 16 lattice. From [29]. 



mass together with its ratio are shown as a function of j3. Assuming a nucleon mass of about 
1 GeV, the diquark mass changes from ~ 500 MeV at /3 = 0.90 to ~ 200 MeV at j8 = 1.10. 
On the small lattice the scale is set by the ground state diquark mass 5(j8) = 7ndiq U ark(j3). The 
phase diagram at zero temperature is then given as a function of the dimensionless parameter jx = 
jj. I mdi quark and the dimensionless lattice spacing a. In Fig. ^ the quark number density is shown. 




Figure 9: Quark number density as a function of the lattice spacing a (left panel) and at j8 = 1 .05 (right 
panel) on a 8 3 x 16 lattice. From [29]. 



Independent of the lattice spacing the quark number density takes it maximum value of n q ^ t = 
2 ■ N c ■ Nf = 14 at large p.. With decreasing lattice spacing, the saturation shifts to larger values of 
chemical potential, indicating that this saturation is only a lattice artifact. The Polyakov loop and 
the (renormalized) chiral condensate show almost the same behaviour as on the larger lattices, see 
Fig. [K]. Furthermore, the onset transition from the vacuum to nuclear matter is studied in Fig. 11. 
At /io ~ 0.5, a transition in the quark number density (left panel) is observed. The value of the 
onset does almost not depend on the lattice spacing, indicating that at smaller values of p. finite 
size effects are less important than for larger values of the chemical potential. In the right panel, 
the transition (shaded region) is compared to half of the diquark mass, and a clear coincidence is 
visible. This indeed verifies that G2 QCD possesses, as advertised above, the silver blaze property 
[43] for baryon chemical potential, i.e. half of the mass of the lightest bound state carrying baryon 



10 



G2 gauge theories 



Bjorn H. Wellegehausen 




Figure 10: Polyakov loop (left panel) and chiral condensate (right panel) on a 8 3 x 16 lattice as a function 
of chemical potential jl and lattice spacing a. From [29]. 




number is a lower bound for the onset transition to nuclear matter. With decreasing lattice spacing 
a, a plateau develops for jSn(<5) < p. < p-i(a), where the quark number density is almost constant. 
For jl > /ii (a) it starts again to increase until it saturates at /i sat . 

6. Conclusions 

Concluding, G2 QCD is a highly interesting arena to investigate both conceptual and practical 
questions. Conceptually, it has already taught us that the center of the gauge group is far less 
relevant than originally anticipated. Most of the salient features of Yang-Mills theory are also 
present for this case with trivial center. It can thus be expected that many other questions may be 
little affected by the center as well. However, it also taught us that the group structure and matter 
representation is important for the physics. 

Investigating practical applications, which particularly involve benchmarks for models and 
continuum methods at finite densities and low temperatures, is only a newly emerging field. It has 
been shown that this is possible. It was furthermore already found that the low and intermediate 
density regime are quite different from simple systems, confirming the situation in QC2D. To fully 
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control this domain, so important for compact stellar objects, much progress will be needed. G2 
QCD will, almost certainly, play an important role in the support of this enterprise in the years to 
come. 
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